Efficient algorithm for a quantum analogue of 2-SAT 



Sergey Bravyi 
IBM Watson Research Center, Yorktown Heights, NY 10598 
(Dated: February 1, 2008) 



Complexity of a quantum analogue of the satisfiability problem is studied. Quantum fc-SAT is 
a problem of verifying whether there exists a state |^) of n qubits such that its A;-qubit reduced 
density matrices have support on prescribed subspaces. We present a classical algorithm solving 
quantum 2-SAT in a polynomial time, ft generalizes the well-known algorithm for the classical 
2-SAT. Besides, we show that for any k > 4 quantum fc-SAT is complete in the complexity class 
QMA with one-sided error. 



Quantum analogues of classical complexity classes have been studied extensively during the last several years. From 
the practical perspective the most interesting of them is BQP — a class of problems that can be solved on a quantum 
computer in a polynomial time with bounded error probability. Another relevant class is QMA — a quantum analogue 
of NP. By definition, a language L belongs to QMA if a membership x € L can be efficiently verified on a quantum 
computer having access to a witness quantum state \ip x ) playing the same role as a witness bit string in the definition 
of NP, see pj. It is widely believed that QMA is strictly larger than NP. 

There are only a few natural problems known to be QMA-complete. The results of the present paper concern the 
problem 2-local Hamiltonian introduced by Kitaev, see 0, • It is a suitably formalized version of a fundamental 
problem in the quantum many-body physics: evaluate the ground state energy of n-qubit Hamiltonian which can be 
represented a sum of two-qubit interactions, H = — ^2 a< bHa,b- 



A theorem proven by Kitaev, Kempe, and Regev at asserts that the problem 2-local Hamiltonian is QMA- 
complete. One way to interpret this result is to identify a classical analogue of the 2-local Hamiltonian problem. If 
one assumes that all interactions H a ,b are diagonal in the standard |0), |1) basis of n qubits and that their eigenvalues 
may be either or 1, the 2-local Hamiltonian problem becomes equivalent to MAX-2-SAT 2]. Indeed, in this case a 
ground state of H can be chosen as a basis vector \x) corresponding to some n-bit string x. An eigenvalue of -ff a ,b 
on \x) is a Boolean function of bits x a ,Xb that can be represented as a conjuction of clauses x a V Xb, {^x a ) V Xb, 
x a V (^Xb), (^x a ) V {~^Xb) (non-trivial interactions involve one, two, or three clauses). Accordingly, one can replace H 
by a list of 2-bit clauses, while the ground state energy of H is equal to the maximum number of clauses that can be 
satisfied simultaneously. 

This quantum-to-classical mapping allows one to attribute QMA-completeness of the 2-local Hamiltonian problem 
to NP-completeness of MAX-2-SAT. On the other hand, classical problem 2-SAT (verify that all 2-bit clauses in a 
given list can be satisfied simultaneously) can be solved in a linear time see also Appendix B. It rises the following 
natural questions: 

• What is a quantum analogue of 2-SAT ? 

• Does there exist a poly-time algorithm (classical or quantum) solving quantum 2-SAT ? 

To define quantum 2-SAT we shall substitute a 2-bit clause C a ^b involving some particular pair of bits x a ,Xb (for 
example, C a ,b =i a V (^Xb)) by a pure two qubit state \4> a ,b) £ C 2 (g) C 2 . Quantum analogue of a property "n-bit string 
x satisfies the clause C 0i b" is a property "quantum state |^) of n qubits has reduced density matrix p a ^ orthogonal to 
|0a,b)"- The latter means precisely that {4>a,b\PaA < t ) a,b) = 0, where p a ,b = Tr c ^ Qi b |\&)(^|. Such a generalization seems 
justified since a two-bit clause C a ^b is not satisfied only for one configuration of variables x a , Xb (in the example above 
x a = and Xb = 1), so a satisfying assignment x must be orthogonal to this configuration. In general, there might 
be several clauses involving a pair of bits a, b. The corresponding two-qubit states \4> a ,b) span a forbidden subspace 
M a ,b C C 2 <g> C 2 , and a satisfying assignment |*) must satisfy (<£|p ,&|<£) = for any \<j>) e M a ,&• Denoting U a: b 
an orthogonal projector onto M a ,b multiplied by the identity operator on all other qubits, we obtain an equivalent 
condition n aj f, = 0. Thus quantum 2-SAT can be defined as follows: 

Input: An integer n and a family of 1- qubit 'projectors {n a j,} ; 1 < a < b < n. 

Problem: Decide whether there exists a state \^) of n qubits such that Tl a ,b \^) = for all 1 < a < b <n. 

Note that there is no loss of generality in imposing constraints for every pair of qubits, since some of the projectors 
Tl a ,b may be equal to zero. 



I. INTRODUCTION 
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The main conclusion drawn in the present paper is that quantum 2-SAT is not harder than its classical counterpart. 
An efficient algorithm for quantum 2-SAT is described in Section^ The key ingredient of the algorithm is Lemma"]] 
that allows one to generate new 2-qubit constraints from existing ones according to simple local rules. It can be 
illustrated by the following simple example. Suppose a state |\&) satisfies constraints IT, 2 \*&) = 112,3 |^) = 0, where 
III, 2 = I^X* - ! 1,2^3, and n 2 , 3 = /i®|* _ >(* _ |2,3 are projectors onto the singlet state, |* _ > = 2" 1 / 2 (|0, 1) — 1 1, 0) ) . 
Taking into account that / — IT. 2 and / — 112,3 are projectors onto the symmetric subspace of qubits 1, 2 and 2, 3 
respectively, we conclude that \^} is invariant under swaps of qubits 1 <-> 2 and 2 <-> 3. Thus is also invariant 
under a swap 1 <-> 3. Therefore we can generate a new constraint IT 3 \^f) = 0, where IT 3 = | \T/ ) 1 1 ,3 <8> -fa- 
Lemma [D generalizes this observation to projectors onto arbitrary two-qubit states (not necessarily entangled). 

The algorithm is defined inductively, such that it allows one to reduce ?i-qubit 2-SAT to (n — l)-qubit 2-SAT in 
a polynomial time. Each step of the induction proceeds as follows (for a detailed analysis see Section [nj). For each 
pair 1 < a< 6 < n we look at the rank r a ,b of the projector II a ,6 (the dimension of the forbidden subspace M. a ,b)- 
If i" a ,b = 4, the problem has no satisfying assignments. If r a> b = 3, a state of qubits a, b is completely specified by a 
constraint II a ,fc |\I>) = alone, so one can exclude them from consideration reducing n to n — 2. If r a .b = 2, qubits 
a and b can be merged into a single logical qubit supported on the zero subspace of II a ,{,. It reduces n to n — 1. If 
r a ,b < 1 for all a and b, we invoke Lemma"!] to generate new constraints. The process of generating new constraints 
terminates whenever r a ^ > 2 is encountered at some pair of qubits (which reduces the total number of qubits), or 
if we arrive to what we call a complete set of constraints. We prove then that any instance with a complete set of 
constraints has a satisfying assignment (moreover, it can be chosen as a product of one-qubit states). The running 
time of the whole algorithm is 0(n 4 ). 

A satisfiability problem in which clauses involve fc bits is known as fc-SAT. Quantum analogue of fc-SAT can be 
specified by a family of fc-qubit projectors {lis}, where S C {1, ...,«} runs over all subsets of cardinality fc. One 
needs to verify whether there exists an n-qubit state |\E') such that lis |^) = for all S. Since fc-SAT is known to 
be NP-complete for fc > 3, one can ask whether quantum fc-SAT is QMA-complete for fc > 3? Before addressing this 
question it should be noted that the definition of quantum fc-SAT given above is not quite satisfactory. Indeed, there 
might exist instances which have no exact satisfying assignments, meanwhile having an approximate solutions |\&) such 
that equations lis |\&) = are satisfied with an exponentially small error. To exclude such cases from consideration 
we shall introduce a precision parameter e, e > l/n Q , a — 0(1), separating positive and negative instances. The 
definition of quantum fc-SAT given above must be complimented as follows: 

Input: An integer n, a real number e — fl(l/n a ), and a family of k-qubit projectors {lis}, S C {1, . . . , n}, \S\ = fc. 
Promise: Either there exists n-qubit state |\&) such that lis = for all S, or ^g^^sV^) > e for all 
Problem: Decide which one is the case. 



To analyze complexity of quantum fc-SAT we introduce a class QMAi. It is defined in the same way as QMA with 
the only difference that for positive instances the verifying circuit accepts with probability one, see Section 11*111 for 
more details. The definition implies that QMAi C QMA. 

We prove that quantum fc-SAT with the promise as above is QMAi-complete for fc > 4. This result is quite 
unexpected since in general a subspace spanned by satisfying assignments of quantum fc-SAT lacks a description in 
terms of Z-qubit projectors for any I < fc. So one could expect the complexity of quantum fc-SAT to grow with fc. Our 
result shows that this is not the case, at least for fc > 4. Whether or not quantum 3-SAT is QMAi-complete remains 
an open question. 

The fact that quantum 5-SAT is QMAi-complete easily follows from a possibility to represent the "computational 
history state" associated with a quantum circuit as a common zero vector of 5-qubit projectors, see 0. A contribution 
reported in the present paper is essentially a reduction from fc = 5 to fc = 4. A mapping from a quantum circuit 
to a family of 4-qubit constraints that we use is a significantly simplified version of a construction proposed in 
(Section 4) for adiabatic quantum computation. 

Remark: Input data for all problems discussed in the paper involve linear operators with complex matrix elements. 
To deal with exact equalities one has to use an appropriate exact representation of complex numbers. A good choice is 
algebraic numbers of bounded degree over the field of rational numbers (roots of polynomials with rational coefficients). 
All common linear algebra tasks for operators whose matrix elements are algebraic numbers can be solved efficiently, 
see books 0,0 for the subject. 

The rest of the paper is organized as follows. Efficient algorithm for quantum 2-SAT is presented in Section [n] 
(for the sake of completeness we outline the standard algorithm solving classical 2-SAT in Appendix B). QMAi- 
completeness of quantum fc-SAT, fc > 4, is proved in Section IlIII A technical lemma needed for this proof concerning 
universality of three-qubit quantum gates with matrix elements from a fixed field is placed in Appendix A. 
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II. EFFICIENT ALGORITHM FOR QUANTUM 2-SAT 



Let {n a ,(,}, 1 < a < b < n be an instance of quantum 2-SAT defined on n qubits. Without loss of generality n > 3. 
Obviously, these data can be encoded by a binary string x of length 0(n 2 ). We shall construct an algorithm that 
takes x as input and outputs one of the following 

• Output 1: x has no satisfying assignments. 

• Output 2: A list of one-qubit states {l^j}} such that \4>i) <£> ■ ■ ■ <8> \ip n ) is a satisfying assignment for x. 

• Output 3: An instance y of quantum 2-SAT defined on n — 1 or smaller number of qubits, such that y is 
equivalent to x. 

The algorithm runs in a time 0(n 3 ). By applying it inductively 0(n) times one can reduce x to an equivalent instance 
involving a constant (say, n = 2) number of qubits which can be solved directly. It should be noted that although the 
algorithm allows one to construct one particular satisfying assignment (SA) for positive instances, this SA might not 
be a product of one-qubit states (the reduction corresponding to Output 3 may transform a product SA for y into 
entangled SA for x). Denote r a b a rank of the projector H a ,b- Without loss of generality < r a & < 3. 

Definition 1 An instance of quantum 2-SAT is called homogeneous iff ' r a ,b < 1 for all a,b. 

(Note that if r a .b — 0, there is no any constraint for a pair of qubits a, b.) The first step of the algorithm is to verify 
whether x is homogeneous. If it is not, x can be transformed into an equivalent instance y on a smaller number of 
qubits (producing Output 3) as explained in Subsections A,B. If x is homogeneous (which is the most interesting 
case) one has to proceed to Subsection C. 



A. Projectors of rank three 



Suppose that r a ^ = 3 for some pair of qubits a, 6, i.e., 

n , 6 = i-\<P){ct>\, |^)ec 2 ®c 2 , {<p\<p) = i. 

An eigenvalue equation Yl a ,b I*) = implies that any SA has a product form \ty) = \4>) a b ® |*')> where I*') is some 
state of the remaining n — 2 qubits. Let us show that |^) is a SA iff |^') is a SA for a new instance y of quantum 
2-SAT defined onn-2 qubits. Indeed, for any projectors P, Q one has 

PI*) = |*), r p|*> = |*), 

0|*) = |*), ltt \PQP\y) - |*). [L) 
Let us choose any pair of qubits c, d and set 

p = i-n a , b = \4>)(cf>\ab, Q = i-n c . d . 

If the pairs (a, b) and (c, d) do not overlap, the right-hand side of Eq. (|TJ is equivalent to 

|*> = |0)a 6 ®|*') 1 n c , d |*') = o. 

On the other hand, if the pairs (a, b) and (c, d) overlap at one of the qubits, say a — c, then the right-hand side of 
Eq. 0} is equivalent to 

\*) = m ab ®m, Q d |f') = i*'), 

where Q is a one-qubit self-adjoint operator implicitly defined by 

\4>) (c^lab {I - II M ) \4>) {(f>\ ab = |0) (4>\ab <8 Qd- 

Thus any 2-qubit constraint imposed on |*) is equivalent to a 1-qubit or 2-qubit constraint imposed on I*'). In other 
words, I*) is a SA for x iff |*') is a SA for a new instance of quantum 2-SAT y defined onn-2 qubits. 
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B. Projectors of rank two 



Suppose that r a f, = 2 for some pair of qubits a, b. We shall argue that qubits a, b can be merged into a single logical 
qubit. It will yield a new instance y of quantum 2-SAT defined onti-1 qubits equivalent to x. 
Indeed, let V : C 2 — > C 2 ® C 2 be an isometry such that 

IIa,6=/-(^ t )«6, V f V = I. 

The constraint n aj h |\E f ) = is equivalent to saying that 

I*) =V C \*'), 

where is some (n — l)-qubit state and the label c refers to a qubit resulting from merging of a and b. The fact 
that V is an isometry implies that for any pair of qubits /, g 

n /)9 |*) = o iff v c tn /lfl v c |* / ) = o. (2) 

If the pairs (a, b) and (/, g) do not overlap, the right-hand side of Eq. J5J is equivalent to il/ ig \^f') = 0. On the 
other hand, if the pairs (a, b) and (/, g) overlap at one of the qubits, say f = a, then the right-hand side of Eq. J2J is 
equivalent to 

Q c , ff |*')=0, where Q c<g = Tl a , g V c . (3) 

Obviously, Q c g is a self-adjoint operator that acts only on two qubits c and g. 

Thus any 2-qubit constraint imposed on is equivalent to a 2-qubit constraint imposed on \^'). In other words, 
is a SA for x iff \W) is a SA for a new instance of quantum 2-SAT y defined on n — 1 qubits. 



C. Homogeneous 2-SAT 

From now on we can assume that x is a homogeneous instance, i.e., all non-zero projectors have rank one. It will 
be convenient to describe a configuration of non-zero projectors by a graph G — (V,E), such that V = {1, . . . ,n}, 
and (a, b) £ E iff Ila.;, ^ 0. Now we can reformulate quantum 2-SAT as follows. 

Input: A graph G = (V, E) with n vertices and a list of rank two tensors </>(° ,b ' assigned to edges (a, b) e E. 
Problem: Decide whether there exists a non-zero tensor ip of rank n such that 

E & b ^ u ...,a n = for any (a,b) € E. (4) 

Comments: One should not confuse a rank of a tensor (number of its indexes) with a rank of a projector discussed 
in previous subsections. From now on we shall omit the summation symbol assuming that all repeated indexes are 
contracted. Indexes of all tensors take values and 1. To avoid clutter in formulas, sometimes we shall identify 
tensors of rank two with 2x2 matrices and use matrix/vector multiplication. Let us also agree that 4>^^ — <f>^'a ■ 
We shall use a symbol e for the fully antisymmetric tensor of rank two, i.e., 




The following important observation provides a simple local rule for generating new constraints. 
Lemma 1 Let <p, 9 be arbitrary tensors of rank two and tp be a tensor of rank three such that 

<t>a,0 ^a,0^ = and 8f3 n ij) a ,p,i = 0. 

Then ip also obeys 



u a ,~f ipa,0,-y = 0, where to an = <f> at p ep t s 9s tl . 
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Using matrix representation of rank two tensors, one has simply oj — (be 8. 
Proof: Indeed, 

Waj VV0,7 = ^Q,0,7 {4>afi #1,7 — 0q, 1^0,7) 

= — <t>n,l (^a,l,7 ^1,7 + "0q,O,7 #0,7) = 0- 

Analogously one proves that u) a ,-y ^0,1,7 = 0. 

□ 

A straightforward generalization of the lemma is this. Suppose ib is a tensor of rank n such that 

^l^ 1 ,..,an=0 and ^L^W--^=° ( 5 ) 
for some tensors <b^ a,b ^ and 0( b,c ) of rank two and some integers a ^= b =/= c. Then ib also obeys a constraint 

w&l i> ai ,..., an = 0, (6) 

where 

Consider an instance a; of 2-SAT specified by a graph G = (V, £*) and tensors <j)( a ' b \ (a, b) 6 £\ We can try to use 
Lemma ^and its corollary Eqs. i|tjl7fl to generate new constraints on ip from the existing ones. Indeed, consider any 
pair of edges (a, b) £ E and (b,c) 6 _E. Let w( a ' c ) be a tensor defined in Eq. Q. One has to consider the following 
possibilities: 

1. uj( a < c ) = 0. We get no additional constraints on tp. 

2. o;( a ' c ) ^ 0, but (a, c) 6 i? and cj( q ' c ) is proportional to <p tya,c \ Again, we get no additional constraints on ib. 

3. a/ a ' c ) ^ and (a,c) <^ E. In this case Eq. © provides a new constraint on tp. Let us add (a,c) into the list of 
edges E and set </>( a ' c ) = uj( a - c ) . 

4. a/ a ' c ) 7^ 0, (a, c) 6 E 1 , and cj( a ' c ) is not proportional to cb^ a - c \ In this case we get two independent constraints 
for the pair of qubits a and c. They leave us with a two-dimensional forbidden subspace for this pair of qubits. 
Thus we can apply ideas of the previous subsection to merge a and c into a single logical qubit and obtain a 
new instance of quantum 2-SAT with n — 1 qubits which is equivalent to the original one. 

Let us keep trying to generate new constraints by probing different pairs of edges until we either encounter the 
case 4 reducing the number of qubits from n to n — 1 (the algorithm terminates with Output 3), or after 0(n 3 ) steps 
we arrive to a homogeneous instance with a complete set of constraints which is defined below. 

Definition 2 A graph G — (V,E) and a family of rank two tensors (b^ a ' h \ (a,b) S E constitute a complete set of 
constraints iff for any pair of edges (a,b) £ E and (b,c) £ E one of the following is true 

• If (a, c) £ E then cb^ e (b^ is either proportional to </A a,c ' or zero, 

• If (a, c) i E then (b^ e (b^ = 0. 

This definition just says that any attempt to generate a new constraint on ip using Lemma ^would fail. 

Lemma 2 Any homogeneous instance of quantum 2-SAT with a complete set of constraints has a satisfying assign- 
ment. It can be chosen as a product of one-qubit state. 

Remark: A proof given below is constructive. It allows one to find a product satisfying assignment in a time 0{n), 
n = \V\. 

Proof: Let {4>( a ' b >}, (a, b) £ E be a family of rank two tensors constituting a complete set of constraints. We shall 
construct a satisfying assignment ip which is a product of tensors of rank one: 

ib - ib {1) ib {2) • • • ib {n) 

Vai,Q2,...,a rl — V ai Va 2 Tan ' 
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Let us start from choosing an arbitrary i/j^K Denote nbgh(a) C V a set of all nearest neighbors of a vertex a. The 
next step is to assign states to all vertices a £ nbgh(l) according to = e (0( 1 >°))- r tp^- 1 ', or in tensor notations, 

i> { : ] =ea,p<$fi> [ J\ aenbgh(l). (8) 
Let us verify that all constraints on edges incident to the vertex 1 are satisfied. Indeed, 

because e is an antisymmetric tensor. This is equivalent to the desired constraints 

<t>aul a ^i,..,a, = for any a G nbgh(l) . 

The next step is to verify that all constraints on edges connecting two vertices from nbgh(l) are automatically satisfied. 
Indeed consider any pair of edges (1, a) € E and (l,c) £ E such that (a, c) 6 £ as well. Then by definition of a 
complete set of constraints we have (j)^ a ^ = (fM'^> e<jy- x,c > (up to some overall factor). Since we have already fulfilled 
the constraints on the edges (a, 1) and (1, c), we can apply Lemma^(with b = 1) to infer that 

0i a V c) V>i a) ^ c) =o. 

This is equivalent to the desired constraints 

<t>a^L ^<w,...,a„ = for an y a ' c e nbgM 1 )- ( a - c ) € 

Let us split the set of all vertices V into two subsets, V — V c iose U Vf ar , where V c iose consists of the vertex 1 and its 
nearest neighbors, while Vf ar consists of all other vertices (those having distance two or greater from the vertex 1). 
All vertices in V c i ose have been already assigned a state. This assignment satisfies the constraints on all edges having 
both ends in V c i ose . We shall now prove that this assignment also satisfies the constraints on any edge having one end 
in V c i ose and another end in Vt ar , for any choice of assignments in Vf ar . 

Indeed, consider any edge (b,c) £ E, such that b £ nbgh(l) and c £ Vfar- Since (l,c) ^ E, by definition of a 
complete set of constraints we have 

0(1,6) e ^(i.c) =a ( 9 ) 
Let us verify that the constraint on the edge (6, c) is fulfilled for all choices of an assignment ip^ c ', i- e -> 

= for all V (c) - (10) 



It suffices to verify that 



Indeed, from Eq. JSJ we have 



$; 7 c) # = o. (ii) 



Therefore 

see Eq. Thus we have proved Eq. Qll|). and therefore Eq. (|l(Jfl . This is equivalent to the desired constraints 

<'l V»ax,...,a» = for any 6 G nbgh(l), c^V (12) 

for any choice of an assignment at the vertices c € V/ or . 

The only constraints that we have not verified yet are those sitting on edges having both ends in Vf ar . One can 
easily check that a family of rank two tensors {</>( a,b )}, a, b £ Vf ar , (a, b) £ E, constitutes a complete set of constraints 
(since its definition is local). Thus we can keep applying the same algorithm to assign states to vertices from Vf ar . 
In this way we will finally end up with a SA for all n qubits. 

□ 
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III. QUANTUM 4-SAT IS QMAi COMPLETE 

Let us start from defining a complexity class QMAi. Throughout this section all quantum circuits are assumed to 
have a data input register R in , a witness input register R w u, and an output register R ou t- We shall assume that all 
qubits of the data input register R in are set initially in |0) state. Final measurement involves measurement of each 
qubit from R out in the standard {|0), |1)} basis. For any quantum circuit U and witness input state \ip w u) define an 
acceptance probability AP(U,ip W i t ) as a probability for the final measurement to yield the outcome in every qubit 
of R out . 

Let T C C be a field that we use for exact representation of complex numbers (see the remark in Introduction). 
We shall consider a gate set Q consisting of all three-qubit unitary operators whose matrix elements in the standard 
basis belong to T (we include three-qubit gates into Q just to simplify the proofs). 

Definition 3 A language L — L yes U L no C S* belongs to QMAi iff there exists a polynomial p and uniform family 
of quantum circuits {U(x)}, x € L, such that U(x) has a length at most p(\x\), U(x) uses only gates from Q, 
\R m \ + \R W it\ + \Rout\ < P(\x\), and 

• If x € Ly es then there exists a witness state \tpwit) such that AP{U{x),^p W it) = 1; 

• If x € L no then for any witness state \ip W it) one has AP{U{x) 1 ip W it) < 1 — l/p(|x|). 

Since there is a polynomial gap between acceptance probabilities corresponding to positive and negative instances, 
we conclude that QMAi ^ QMA. From now on we shall assume that for any instance of quantum fc-SAT projectors 
Us have matrix elements from the field T . 

Lemma 3 Quantum k-SAT belongs to QMAi f 0T ony constant k. 

Proof: Let {lis}, S C {1,2,..., n}, \S\ — k be an instance of quantum fe-SAT. We shall construct a quantum circuit 
that measures eigenvalues of Ilg and writes the outcomes to output register R ou t- More strictly, for each subset £>, 
l^l = fc, introduce an auxiliary qubit q(S) initially prepared in the |0) state, and consider a unitary operator 

u s = n s ® of (s) + (/ - u s ) ® i q(s) . 

It acts on a witness input register R w a of size n and the auxiliary qubit q(S). Note that Us acts non-trivially only 
on a subset of fe + 1 qubits. According to Lemma 03 Appendix A, the operator Us can be exactly represented by a 
quantum circuit of a size poly(2 k ) with gates from the set Q . 
For any input witness state \ip W it} one has 

u s Wwit) ® |o) = n s \ip wit ) ® |i> + (i - n s ) IVw) <8 |o). (13) 

Denote ps a probability to measure the qubit q(S) in the state |1) after the application of Us- From Eq. i(T3|) one can 
easily get 

Ps = (Vwl n s|Vw)- 

Let us define the output register R ou t as a collection of all ancillary qubits q(S), one qubit for each subset 5. Choose 
some order on the set of subsets S and define a verifying circuit as U — Y\ s Us- By construction, U accepts on a 
witness state \ip W it) iff all the qubits q(S) have been measured in |0). A probability of that event can be bounded as 

1 ~J2 Ps <AP{U,^ wit ) < l-maxpg. (14) 

s 

If there exists a satisfying assignment, one can choose \ip w u) such that ps = for all 5 1 , and thus AP{U,ij} W it) = 1. 
On the other hand, if there is no satisfying assignment, we are promised that for any witness state ^2 S PS — e ' wnere 
e = l/n a , a — O(l), see Introduction. Therefore 

max p s > > n' k e. 

Combining it with Eq. (|14|l we get AP(U, i/Wt) — 1 — l/n k+a for any witness state \ip W it)- Now one can easily choose 
a proper polynomial p to meet all the requirements of Definition |21 
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□ 



Lemma 4 Quantum k-SAT is QMAi -complete for any constant k > 4. 

Proof of Lemma [IJ It suffices to prove that quantum 4-SAT is QMAi-complete. Let U = Ul ■ ■ ■ U% U\, Uj € Q, 
be a quantum circuit operating on TV qubits with input data and witness registers Ri n and R w u, and output register 
R ou t- Denote \ipin) = |0, . . . , 0) and \ipwit) the initial states of Ri„ and R w u- Without loss of generality we can assume 
that N = \Rin \ + l-Rtuitl- The collection of N qubits Ri n U R w u there the computation goes will be referred to as a 
computational register. It is described by a Hilbert space 

H = fC 2 ) 0JV 

Let us define L + 1 intermediate computational states 

I Qo) = \ipin) ® |VW and I Qt) = U t \Qt-i), t =l,...,L. 

Introduce an auxiliary Hilbert space 

^clock = (C 4 )® L 

that describes a clock register composed of L four-dimensional clock particles, one particle associated with each gate 
in the circuit. Each clock particle can be represented by a pair of qubits, so if we want to keep all projectors 4-local, 
we can afford only projectors acting on n c clock particles and n q computational qubits, where 

(n c ,n,)€ {(2,0), (1,2), (0,4)}. 

Later on we shall define 2L orthonormal legal clock states 

, \C[), . . . , \Cl)i \C l ) £ Hciock- 

Our first goal is to design a Hamiltonian H(U) such that 

• H(U) acts on a Hilbert space H = H c iock ® H com p, 

• H(U) is a sum of 4-qubit projectors; 

• The zero-subspace of H(U) is spanned by vectors 

L 

m wi t)) = ^\Ct)®\Qt-i) + \C' t )®\Qt), \Qo) = \il)in)®Wwit)- (15) 
t=i 



The state \Q(ip W it)} represents the whole history of the computation starting from a witness state \ipwit)- Note that 
the gate Ut is applied when the clock register goes from the state \Ct) to the state \C' t ). On the other hand, as the 
clock register goes from \C' t ) to \C t +i), the state of the computational register does not change. 

Now let us specify the legal clock states. Basis states of a single clock particle will be denoted as \u) — unborn, \d) 
- dead, \al) — active phase 1, \a2) — active phase 2 (this terminology is partially borrowed from [j]. As the time 
goes forward, each clock particle evolves from the unborn phase to active phase 1, then to active phase 2 and finally 
ends up in the dead phase. Below we list the eight legal clock states for L = 4 as an example. 

\Ci) = \al,u,u,u), |C 3 ) = \d,d,al,u), 

\C[) = \a2,u,u,u), jCg) = \d,d,a2,u), . . 

\C 2 ) = \d,al,u,u), \d) = \d,d,d,al), 

\C 2 ) = \d,a2,u,u), \C' A ) = \d,d,d,a2). 

In general, legal clock states are defined as basis vectors of 7i c iock that obey the following constraints: 

1. The particle 1 is either active or dead, 

2. The particle L is either active or unborn, 

3. There is at most one active particle, 

4. If a particle j is dead for all 1 < k < j a particle k is dead, 
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5. If a particle j is unborn then for all j < k < L a particle k is unborn, 

6. If a particle j is dead then a particle j + 1 is either dead or active. 

One can easily check that there are only 2L basis states satisfying these rules. Each of these states has exactly one 
active particle (in cither active phase 1 or active phase 2) at some position 1 < j < L, dead particles at positions 
k < j and unborn particles at positions k > j, 

\Cj) = | d, . . . , d, ol, u, . . . , u), \CA — | d, . . . , d, a2, u,...,u). 



3 i _ ____ _____ 

Denote 

W lcgal = Linear Span(|C_), \C[), . . . , \C L ), \C' L )) C W clock 
a subspace spanned by the legal clock states. Let us define a Hamiltonian iJ c iock acting on the space Wdock as 

6 

u \ " itU) 

-"clock — 2_^i -"clock' 

^ciock = \d)(d\ L , 

^ciock = E (Ial><al| + |a2)<a2|) j ®(|al)<al| + |a2)<a2|) fc , 

l<j<k<L 

H^L = E (\al)(al\ + \a2)(a2\ + \u)(u\)^\d)(d\ k , 

l<j<k<L 

^ciock = E l«>H-®(|al><al| + |a2)<a2| + |d><d|) fc) 

l<j<k<L 

^ciock = E \d)(d\ 3 ®\u)(u\ 3+1 . 

1<3<L-1 



The six terms of H c \ 0( ± correspond to the six constraints listed above, so that 

IV>)eWiegai iff #clock !?/>)= o. (17) 
Note that .f/ciock is a sum of 4-qubit projectors. Denote II an orthogonal projector onto a subspace Hi ega ,i (g> H C om P , 

L 

(\n.\m.\ _i_ m'\/n'i\ ^ , 

f comp- 



n = J2(\c j )(c j \ + \c' j )(C j \)®i cl 



Define a Hamiltonian iJj n it as 



ffinit = |al)<al|_® ( E |l><l|b ) • ( 18 ) 

The only legal clock state having \al) in the first position is \C\). Thus i?i n it penalizes any qubit of the input data 
register R in for being in a state |1) provided that the clock register's state is |Ci). A restriction of i?i n i t onto the 
subspace Hi ega i ® H CO mp is therefore 

nfr init n = |Ci)<c7i|® ( E |i><i|b) - (19) 

\beR in / 

Note that H m i t is a sum of 3-qubit projectors. 
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Define a Hamiltonian H pmp as 



L-l 



where 



and 



-^prop,t — 2 



Hpmp — -ffprop,* + H prop t , (20) 
t=i t=i 

(|ol)(ol| + |a2)(a2|) t (g)/ com p- \a2){al\ t ®U t - \al) (a2\ t ®u(\ , (21) 



#prop,t = g (|o2,u)<a2,u| + |d,ol)(d,ol| - |d,ol)(o2,«| - |a2, al|) M+1 ® / comp . (22) 

The operators H pmPtt and H' pmp t are 4-qubit projectors. By inspecting the example Eq. Ijl6(l one can easily check 
that 

• The only legal clock state having al at position t is |C t ); 

• The only legal clock state having a2 at position t is \C' t ); 

• The only legal clock state having (d, al) at positions t, t + 1 is |Ct+i); 

• The only legal clock state having (a2, it) at positions t, t + 1 is |C t '). 
Thus a restriction of H proPtt onto the subspace Tiiegai ® H CO mp is 

1 
2 

Analogously, a restriction of H' pt onto the subspace Ti.i ega i ® H com p is 



(|C t )<C t | + |CJ)(C?|) ® I comp - |C t ')(C t | ® 17*- |C t )(C t '| ® t7?l . (23) 



n#;, oPit n = \ (\c' t )(C t \ + \c t+1 ){c t+1 \ - \c t+1 ){C t \ - |c t ')(c t+1 |) ® / comp . (24) 

Define a Hamiltonian 

-ff(f) = ffinit + ffciock + ^prop- 
Combining Eqs. (|19I23I24|) one can easily check that the zero-subspace of H(U) is indeed spanned by computational 
history states \Sl(tp W it)) ■ 

One remains to introduce an extra term into H(U) that is responsible for the final measurement of R ou t- Define 



\a2){a2\ L ®( J2 I 1 ) (Mb) 



The only legal clock state having a2 at the position L is \C' L ). Thus a restriction of H out onto the subspace Hi ega \ 

icomp 



E WW 



ni/ out n = |ai)(c?i|® 

Therefore H out penalizes any qubit of the output register for being in the state |1) provided that the clock register is 
in the state \C' L ). Besides, H oni is a sum of 3-qubit projectors. 
We summarize that a Hamiltonian 

H = H init + i?dock + -ffprop + H out (25) 

has zero ground state energy iff there exists input witness state \ip W it) such that AP(U,i/) W it) = 1. By construction, 
H operates on N + 2L qubits, and it can be represented as a sum of 4-qubit projectors: 

H = Y,n s , S C{1,...,N + 2L}, \S\=4. (26) 
s 

Thus the property "H has zero ground state energy" is equivalent to the quantum 4-SAT {n^} having a satisfying 
assignment. 

Let L = L yes U L no be a language from QMAi, x € L be a binary string, and U{x) be a verifying circuit see 
Definition [3J Using the majority voting to amplify the gap in acceptance probabilities, see Q , we can assume that 
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• If x € L yes then AP(U,4> w n) = 1 for some input witness state \tp w u}, 

• If x e L no then AP(U,ip wit ) < e, e = l/p(|a;|), for all |Vw), 

where U is a circuit implementing several copies of U(x) and the majority voting (obviously it can be realized using 
the gate set Q). The polynomial p above may be different from the one in Definition [3J 

To complete the proof of the lemma we have to show that for negative instances, x € L no , the ground state energy 
of the Hamiltonian Eq. (|25l) is not too small, i.e., (^ r |i?|^ r ) > l/g(|x|) for any | \E r ) , where q is some polynomial. It can 
be done using ideas from |l|. Indeed, a decomposition 

H = ^lcgal ® 'Wcomp ffi ^lcgal ® ^comp 

is invariant under H. Let H' be a restriction of H onto 7ii eg ai <8> H CO mp- One can easily check that H' is exactly the 
Hamiltonian that one would assign to a quantum circuit U using the construction of pj , where U is obtained from U 
by appending the identity gate to each gate of U (if U has a length L then U has a length 2L). It was shown in pj 
that H 1 has the ground state energy at least c(l — -ye)(2L) _3 , where c is a constant and e is defined above. It may 
be only polynomially small in |x|. 

Let H" be a restriction of H onto HjJ; gal ® H C omp- Since any state from Hjj; gal violates at least one constraint from 
-ffciock, the ground state energy of H" is at least 1. 

□ 

Remark 1: One could also try to consider clock particles with only three states: unborn, active, and dead. It is 
possible to design a proper dynamics of the clock register, such that the computational history state can be specified 
by 3-local projectors 0. Each projector involves a triple of particles with dimensions 3x2x2. 

Remark 2: If one modifies Definition [3] such that the acceptance probability corresponding to negative instances is 
1 — e, where e > may be arbitrarily small, the corresponding class would be hard for the polynomial hierarchy, as 
it would contain the exact counting class C = P, see 0- The same remark concerns quantum fc-SAT, k > 4, without 
the polynomial gap promise. 
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APPENDIX A: QUANTUM CIRCUITS OVER A FIXED FIELD 

Lemma 5 Let V be a unitary operator acting on k qubits. Suppose that matrix elements of V in the standard basis 
belong to a field T C C. Then V can be exactly represented by a quantum circuit of size poly(k) ■ 2 2k with three-qubit 
gates whose matrix elements in the standard basis also belong to the field T . 

Proof: The first step directly follows the proof of universality of two-qubit gates, see @- Let us decompose V as 
V = W1W2 ■ ■ ■ Wl, where each operator Wj is a direct sum of 2x2 unitary on a two-dimensional subspace spanned 
by some pair of basis vectors and the identity operator on a subspace spanned by the remaining 2 k — 2 basis vectors. 
Accordingly, L — 2 k ~ 1 (2 k — 1) is the number of such subspaces. Matrix elements of the operators Wj belong to the 
same field T . Using swaps of qubits and o~ x gates one can transform any operator Wj as above into a controlled 
one-qubit gate with one target qubit and k — 1 control qubits, Wj = Ak-i(Uj). Here Uj is a one-qubit gate with 
matrix elements from T . 

Consider a classical unitary operator f2 that reversibly computes logical AND of k bits x\, ... , x^-x, i.e., 

VL I a) <g> \xi,x 2 , ■ • ■ ,Xk-i) = \a® AND(xi, . . .,x k -i)) <8> \xi,x 2 , ■ ■ .,x k -i}- 

Here \a) refers to an ancillary qubit. Since three-bit classical gates constitute universal basis for classical reversible 
computation, we can implement f2 (probably using additional |0) ancillas) by a circuit of size poly(k) with three-qubit 
gates whose matrix elements are only and 1. Now one can implement A/-_x(f7,*) as follows (we label the ancillary 
qubit by A and the target qubit by T): 

(1) Set A to |0), 

(2) Apply VL to A and the k — 1 control qubits, 

(3) Apply A(Uj) to A and T such that A is the control qubit, 

(4) Apply SI -1 to A and the k — 1 control qubits. 

One can easily check that the ancillary qubit ends up in the |0) state, while the k — 1 control qubits and the target 
qubit are acted on by Ak—i(Uj). By composing these circuits for each j one gets a circuit representing V. 

□ 



APPENDIX B: EFFICIENT ALGORITHM FOR CLASSICAL 2-SAT 

Let x = (a?i, . . . , x n ) be a binary string. A Boolean function l(x) is called a literal if l(x) = x a or l(x) — (—<x a ) for 
some 1 < a < n. Given m Boolean functions (clauses) Cj(x) — lj(x) V l'j(x), where lj and l'j are literals, a classical 
2-SAT problem is to determine whether there exists a string x such that Cj(x) = 1 for all j = 1, . . . , m. Below we 
briefly describe a well-known algorithm for a classical 2-SAT that runs in a time 0{n + m). 

Without loss of generality l'j ^ ("'lj) (otherwise Cj(x) = 1 for all x). Consider a directed graph G = (V,E), whose 
vertices are literals, i.e., 

V \X\ j • • • ; X n , X~y , . . . ,X n }, 

and whose edges are pairs of literals (I, V) that appear in the same clause as shown below: 

E = {(1, 1') : Cj = lV H') for some j}. 

Each clause Cj contributes to one or two edges of G. Let us split V into strongly connected components (SCC). (By 
definition, / and /' belong to the same SCC iff there exists a path from I to I' and a path from V to I.) Obviously, one 
can identify all SCCs of G in a time 0{n + m). 

Lemma 6 A satisfying assignment x exists iff for any I £ V vertices I and -^l belong to different SCCs. 
Proof: 

(a) Suppose for some 1,1' £ V there exists a path from I to If x is a SA then l(x) = implies l'(x) — 0. If / and ~^l 
belong to the same SCC, then l(x) = implies -<l(x) = and -<l(x) — implies l(x) — 0. This is a contradiction. 

(b) Suppose that / and ~^l belong to different SCCs for all / £ V. Let us find a SA. Consider a directed graph 
G = (V, E) whose vertices are SCCs of G and whose edges are pairs (S, S') of SCCs of G such that there exists a 
path from S to S'. By definition, G is acyclic. Using the topological sorting algorithm one can order vertices of G 
such that (S, S') € E implies S < S' . It can be done in a time 0{n + m). 
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By definition of G, for any vertices I, I' € V a path from I to /' and a path from —V to ->l exist or do not exist 
simultaneously. Thus for any vertex S of G there exists a unique vertex ->S of G such that ^5 = {/ e V : -1/ G S 1 }. 

For all vertices S of G do the following: If 5" < (->S) then set I = 1 for all / e S and set Z = for all / e (-iS 1 ). If 
S > ( _, <S I ) then set Z = for all Z € S and set Z = 1 for all I e ( _, <S I ). We end up with some assignment x. Let us check 
that Cj(x) = 1 for all j. Indeed, otherwise we would have lj(x) = lj(x) = for some j. Let S and S' be SCCs of lj 
and lj. We have S 1 > and S" > -iS". On the other hand, S < ^S' and S' < -*S. This is a contradiction. Thus x 
is a SA. 



□ 



